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There  is  verified  the  hypothesis  about  multifractal  nature  of  radioactive  contami- 
nation due  to  nuclear  explosions  on  Semipalatinsk  test  site  (STS)  in  Kazakhstan. 
The  fields  of  terrestrial  contamination  have  extreme  high  variability  caused  by  a 
number  of  nature  processes.  There  are  used  the  methods  of  multifractal  formalism 
and  computation  topoiogy.  The  results  of  analysis  support  the  existence  of  mul- 
tifractal scaling,  which  is  different  for  natural  and  man-made  isotopes.  It  means 
that  there  are  ” hot  spots”  of  man-made  Cs  isotope  contamination  dangerous  for 
human  health,  which  can’t  be  detected  by  the  traditional  in  Kazakhstan  technics 
of  measurements  from  air.  Moreover,  the  finding  of  self-similarity  could  be  the 
base  of  new  methods  of  diagnostics  of  large  territories. 


1 Introduction 


There  were  some  470  nuclear  explosions  on  STS  in  Kazakhstan,  of  which  90  ex- 
plosions were  in  the  air,  25  on  the  ground  and  355  underground.  As  a result  a lot 
of  territories  were  contaminated  and  it  has  become  a matter  of  serious  for  people 
with  regard  to  its  consequences  for  the  lives  and  health  ^ . We  analyze  the  resulting 
contaminated  areas  by  examining  the  distribution  of  radionuclides.  The  analysis 
of  the  problem  is  not  trivial  and  is  influenced  by  a number  of  issues,  such  as  wide 
range  of  contamination,  extending  from  meters  to  hundreds  of  kilometers;  presence 
of  a mixture  of  radioactive  isotopes,  due  to  numerous  nuclear  explosions  of  vari- 
ous composition  and  power;  highly  variable  rate  of  migration  of  radionuclides  in 
soils;  the  influence  of  climatic  conditions  on  the  localisation  of  contamination  and 
its  transport  ^ . So,  the  radioactive  contamination  of  the  ground  is  a result  from 
nonlinear  interplay  of  geophysical  fields  which  intervene  over  a large  range  of  scale. 
Consequently,  fields  of  contamination  observed  have  an  extreme  spatial  variability, 
frequently  cited  ’’hot  spots”  or  ’’leopard’s  skin” 

Complex  structure  of  the  field  is  the  main  difficulties  of  it’s  diagnostics.  The 
absence  of  the  field  smoothness  restricts  a choice  of  methods  of  measurements  and 
data  processing.  The  maps  of  contamination  are  designed  using  interpolation  of 
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measurements  and,  consequently,  depend  on  a scale  of  averaging.  As  a result, 
sparse  anomalies  can  be  missed. 

To  detect  all  fluctuations  of  the  contamination  field  measurements  must  cover 
a territory  by  a network  of  overlap  neighborhoods^.  But  it  is  impossible  for  huge 
areas  of  former  nuclear  weapon  testing  site,  spanning  an  area  of  around  ISSOOfcm^ 
and  regions  bordering  this  test  site  which  we  deal  with.  Nowadays,  only  1%  of 
contaminated  territory  has  been  investigated  by  accepted  in  Kazakhstan  methods 
of  measurements  from  air  and  it  is  definitely  not  enough  to  describe  the  radioactive 
contamination  and  to  implement  effective  decontamination  methods.  The  scale- 
symmetry  assumption  is  the  simplest  and  also  only  one  acceptable  in  the  absence 
of  knowledge  in  specific  mechanisms  producing  these  fields.  The  idea  of  scaling 
invariance  can  be  useful  in  practice  and  could  become  a base  for  the  analysis  and 
diagnostics  of  the  areas. 

The  present  study  deals  with  the  area-analysis  of  the  distribution  of  radionu- 
clides in  Kazakhstan.  The  data  is  obtained  primarily  from  former  STS  and,  sec- 
ondly, from  regions  bordering  this  test  site.  Experimental  data  consist  of  radionu- 
clide measurements  made  from  the  air,  on  a grid  of  parallel  lines  covering  a part 
of  territory.  Two  regions,  Karaganda  and  Semipalatinsk,  and  Irtysh  test  site  were 
explored. 

It  is  possible  to  analyze  these  measurements  by  several  complementary  methods, 
among  them  the  fractal  multifractal  and  morphology  methods  Here,  we 
present  multifractal  method  and  method  of  computational  topology 

The  structure  of  this  paper  is  as  follows.  In  Sec.  2 we  recall  some  definitions 
of  multifractal  analysis  and  apply  this  technique  to  the  data  from  flight  paths.  In 
Sec.  3 we  describe  the  method  of  computational  topology  and  present  the  results 
of  applying  this  method  and  multifractal  one  to  data  from  the  Irtysh  test  site.  The 
summary  is  found  in  the  concluding  section. 

2 Multifractal  analysis  of  radioactive  contamination 

We  used  an  approach  in  which  the  contamination  data  along  a path  are  considered 
as  multifractal  random  measures.  Let  us  remind  that  multifractal  spectrum 

of  singularities  a of  Borel  finite  measure  /r  on  a *compact*  set  X is  a function  /(q) 
defined  by  pair  {g,  G).  Here,  g : X [—00,00]  is  a function,  which  determines  the 
level  sets: 


K^  = {x£X  : g{x)  = a}  (1) 

and  produces  a multifractal  decomposition  X : 

X=  [J  (2) 

— oo<oi<oo 

Let  G be  a real  function  which  is  defined  on  Zi  C X such  as  G{Zi)  < G{Z2)  if 
Zi  C Z2.  Then  multifractal  spectrum  is  /(a)  = G{Kg).  Let  g be  determined  as 
pointwise  dimension  of  measure  fi{x)  at  all  points  a;  £ A"  for  which  the  limit 

g = d^{x)=  Urn  log  fi{B{x,r))/ log r, 

r—*0 


(3) 
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exists,  where  r))  is  a ’’mass”  of  measure  in  the  ball  of  radius  r centered  at  x. 

Since,  we  have  chosen  g = can  drop  the  g subscript  from  further  references 

to  K^.  Then  Ka  = {a;  : d^(a:)  = a},  where  exponent  a is  local  density  of  p.  The 
singular  distribution  fi  can  then  be  characterized  by  giving  the  ’’size”  of  sets  Ka 
by  their  Hausdorff  dimension,  i.e.  /(a)  = G{Ka)  = dimH{Ka)-  If  M is  self-similar 
in  some  sense,  /(a)  is  a well  behaved  concave  function  of  a To  estimate  f{a)  for 
path  data  we  applied  the  method  of  the  partition  sum 

Let  X = S = [0, 1]  and  be  an  increasing  sequence  of  positive  integers.  Define: 


i i -h  1" 

A»n 


(4) 


The  component  of  multifractal  decomposition  S is 

log 


Ka  = {x  € \ lim 


logA;„ 


= a}. 


It  is  more  easy  to  analyze  ’’coarse  grained”  Holder  exponents  a 

^i.n  ■ 


log/x(/j,n) 

logfc„ 


Their  variations  are  described  by  Legendre  spectrum 

f{a)  = fi{a)  = inf(ag  - T{q)), 
9 

where 


(5) 

On  of  p on 

(6) 

(7) 

(8) 


and  summation  runs  through  i such  that  7^  0. 

We  analyze  the  radionuclide  contamination  data  measured  from  the  air  for  Kara- 
ganda and  Semipalatinsk  regions.  The  height  of  all  flights  was  about  50  m and  the 
velocity  of  the  plane  was  about  125  km/h.  Spectrometrical  measurements  were  usu- 
ally made  at  time  intervals  of  1 sec.  The  technique  used  is  based  on  measurements 
of  7-quanta  flow  density  of  (1.12  and  1.76  Mev)  to  determine  the  contam- 
ination by  U,  (2.62  Mev)  to  determine  Th  contamination,  (1.46  Mev) 

to  determine  K and  (0.66  Mev)  for  Cs  contamination.  Total  7-activity  is 

measured  in  the  range  of  0.25-3.0  Mev.  The  spectrometer  was  equipped  by  Nal 
element.  The  measurements  have  been  done  at  different  scales,  e.g.,  in  Karaganda 
and  Semipalatinsk  region  the  scale  is  1:1,000,000,  so  the  distance  between  paths 
is  10  kms.  There  are  up  to  11000  data  points  along  the  paths,  with  each  reading 
separated  by  50  meters.  A sample  of  measured  data  can  be  found  in  Fig.  1. 

For  calculation  of  multifractal  spectrum  it  is  necessary  to  transform  contamina- 
tion values  Y into  probabilistic  measures.  As  a rule,  empirical  probability  density 
function  p{Y)  behaves  itself  as  a curve  (1)  shown  in  Fig.  2,  where  function  tails 
fall  down.  The  path  data  Y were  transformed  to  new  variables  Z which  have  rect- 
angular random  probability  density  function  (the  curve  2 in  Fig.  2)  by  solution  of 
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Figure  1.  Data  array  of  Cs,  Semipalatinsk  region  path 


bin 


Figure  2.  The  histograms  of  distribution  contamination  data  (1)  and  probabilistic  measures  (2). 


equation  P{Y)  = P{Z),  where  P is  a probability  function.  In  the  case  of  rectangu- 
lar random  distribution  Z = P{Y).  This  transformation  is  equivalent  to  coordinate 
change  which  saves  a dimension  spectrum  /(a) 

The  range  of  q value  was  chosen  as  [—6,4].  The  Fig.  3 displays  the  set  of 
multifractal  spectra  of  radionuclide  contamination  of  paths.  The  scaling  of  all 
isotopes  obeys  to  multifractal  law.  The  behaviors  of  /(a)  curves  are  different  for 
various  isotopes.  In  Fig.  4 it  is  shown  the  set  of  f{a)  spectra  of  man-made  isotope 
Cs  data  calculated  for  different  paths.  It  can  be  seen  that  scaling  properties  differ 
from  path  to  path,  but  this  scaling  variability  is  not  so  high  as  distinction  of  /(a) 
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a distance  between  adjacent  paths  is  ~ lOfcm  and  effective  strip  of  equipment 
capture  is  300m.  For  the  case  of  measurements  paving  the  territory  we  can 
use  the  method  of  computational  topology  described  in  the  next  section,  which 
complements  multifractal  method. 


3 Computational  topology  of  the  radionuclide  fields: 
disconnectedness  index. 


The  multifractal  formalism  can  be  applied  to  area  contaminations,  where  a support 
of  measure  fJ,{x,y)  is  a compact  set  X = [0, 1]  x [0, 1]  G 3?^.  In  this  case  ^ 

where 


' i 2 + T 

C f 2 + T 

ic  ’ Jc 

U I\,YI  ri'n 

^ k ' k 

(9) 


Unfortunately,  contamination  data  do  not  allow  to  use  square  paving,  as  the  dis- 
tance between  the  data  points  on  adjacent  paths  is  always  larger  than  the  distance 
between  the  neighbouring  points  along  a single  path.  It  could  be  used  the  tessela- 
tion  of  the  surface  by  rectangles  {ci}  In  this  case  probability  measure  on  Ci  is 
averaged  over  rectangle,  such  as,  y,{x,  y)  : x,y  G Ci  presents  itself  step-like  function, 
which  has  a non- zero  constant  value  at  each  rectangle.  The  components  of  ’’coarse- 
graining”  multifractal  decomposition  are  found  by  selecting  of  rectangles  Cj  of 
Holder  exponent  a ± Aa,  where  Aa  is  discrete  step  of  the  measure.  The  coarse 
Holder  exponent  is 


o:(ci) 


logM(ci) 

log  |c| 


(10) 


where  |c|  is  a size  of  Cj.  ” Coarse-graining”  version  f{a)  is 


/g(q) 


lim  lim 
Aa-.0i-.[c| 


logNs{a,Aa) 

log(l/(5) 


(11) 


where  Ns{a,Aa)  denotes  the  number  of  boxes  of  6 size  with  a(cj)  ~ a.  This 
function  can  be  explained  in  statistical  terms  as  probability  of  finding  a(cj)  ~ 
Q and  behaves  roughly  like 


Ns{a)/Ns  = PS  [a{d)  ~ a]  ~ 


(12) 


Here,  Ns  is  the  total  number  of  non-empty  (5-boxes  which  contain  measure,  D 
denotes  the  box  dimension  of  the  support  of  /i. 

For  estimation  of  foi^)  spectrum  it  could  be  used  functional  box-counting 
however  the  covering  by  rectangle  neighborhoods  makes  it  impossible  to  re-scale 
rectangle  edges  uniformly  and  introduces  additional  local  and  global  dimensions 
which  have  no  clear  meaning  Consequently,  we  use  computational  topology  for 
calculation  of  the  box  dimension 

The  component  Ka  of  multifractal  decomposition  is  a subset  {d}  having  the 
same  a±Aa.  The  ’’coarse-graining”  geometry  of  Ka  can  be  described  by  the  rate 
of  growth  a number  of  d £ Ka  with  improving  resolution  and  called  disconnect- 
edness index  7. 


Figure  5.  The  paving  map  of  U contamination  of  Irtysh  test  site. 


A topological  space  X is  connected  if  and  only  if  it  cannot  be  decomposed  into 
the  union  of  two  non-empty,  disjoint,  closed  sets.  If  such  a decomposition  exists 
then  X is  said  to  be  disconnected-  that  is,  if  there  are  two  closed  sets  U and  V 
such  that  U f]V  = 0 and  U^jV  = X.  There  was  reformulated  the  notion  of 
connectedness  in  order  to  make  it  possible  to  implement  a test  for  this  property  on 
a computer.  The  basic  idea  is  to  look  at  the  set  with  a finite  resolution  e and  see 
how  connectedness  changes  as  we  let  e ^ 0.  It  is  said  that  a subset  X of  metric 
space  is  £ — disconnected  if  it  can  be  decomposed  into  sets  that  separated  by  a 
distance  of  at  least  e. 

Let  N^ia)  be  the  number  £-components  by  given  a resolution  £ and  intensity  of 
measure  equals  a ± Aa.  The  disconnectedness  index  7 is  the  following  limit: 


7 = lim  inf 
e-*0 


logN^ja) 
log{l/e)  ■ 


(13) 


The  index  equals  to  Hausdorff  dimension  for  the  simplest  sets,  for  example  for 
the  middle-thirds  Cantor  set.  However,  more  generally  the  disconnectedness  index 
distinguishes  between  sets  with  the  same  box-counting  dimension. 


3. 1 Scaling  properties  of  Irtysh  area  contamination. 

There  have  been  applied  both  multifractal  formalism  and  computational  topology 
to  Th,  K,  U,  Cs  isotope  contamination  of  Irtysh  area  which  is  the  part  of  STS.  Irtysh 
area  was  chosen  for  calibration  of  the  devices  for  spectrometrical  measurements  from 
air  and  was  investigated  more  thoroughly.  The  distance  between  paths  was  50  m, 
effective  strip  of  equipment  capture  was  ~ 50m  and  there  were  1197  measurements 
on  the  ground. 

To  calculate  7 index  a range  of  values  for  each  isotope  have  been  divided  into 
255  levels.  The  contamination  values  r'(cj)  are  expressed  in  the  share  of  the  whole 
range  which  was  reduced  to  1.  In  Fig.  5 a part  of  the  digital  image  of  contamination 
used  for  7 computation  are  shown. 

The  curves  of  disconnectedness  index  (Fig.  7)  of  natural  isotopes  have  plateaux 
with  7 = 2 at  S [0.3  -r-  0,5]  for  Th,  v G [0.3  0.45]  for  K and  v G [0.3  0.35]  for 

U and  the  7 values  decrease  as  v increases.  The  behavior  of  man-made  Cs  isotope 
is  different.  The  index  of  disconnectedness  7 is  equal  to  1.3  at  a small  y and  slowly 
decreases  demonstrating  two  plateaux:  7 = 1 at  t/  G [0.45-^0.7]  and  plateau  7 ss  0.6 
at  y > 0.8  (438  milliCurie/Zcm^).  The  latter  means  the  possibility  of  existence  of 
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’’hot  spots”  of  Cs  contamination  which  lie  in  the  area  of  highest  concentrations 
(>  4:38milliCurie/km?)  and  small  fractal  dimensions  7 < 0.6.  In  general,  7(1^)  can 
be  interpreted  as  ” pre- multifractal”  spectrum  /g(q:).  In  this  case,  the  probability 
of  existence  of  ’’hot  spots”  decreases  with  spot  size  6 as  (x  (cf.  Eq.  12).  It 

is  necessary  to  notice  the  similar  behavior  of  Cs  cumulative  soil  deposition  data  in 
Chechslovakia  after  Chernobyl  accident 

Statistical  properties  of  multifractal  fields  are  a function  both  of  a scale  and 
a dimension  of  a support  of  measurements.  It  is  known,  that  measurement  net 
with  empirical  box-dimension  Dnet  cannot  detect  sparsely  distributed  phenomena 
with  box-dimension  Da  < 2 — Dnet  This  conclusion  follows  from  the  common 
ideas  about  tranversality  It  is  necessary  remember,  that  two  sets  M and  N of 
dimension  Dm  and  Dfj  intersect  transversally  in  Sft’”  if 

codim{M  D TV)  = codimDM  + codirnDM-  (14) 

So,  the  net  with  dimension  Dnet  can  detect  a set  of  dimension  Da  only  if  Dnet+Da  > 
m.  As  a rule,  box-dimension  of  the  net  of  spectrometrical  measurements  from  air  is 
Dnet  'C  2,  consequently,  the  number  of  anomalies  on  contaminated  territories  can 
be  missed. 

In  addition  to  computation  topology  method  there  was  used  multifractal  formal- 
ism for  estimation  of  /(a)  spectra  of  Irtysh  test  site  data.  There  was  constructed 
one  dimension  array  by  linking  all  paths  in  one.  Adjusted  paths  were  ’’glued”  to- 
gether by  taking  into  account  their  orientation  and  adding  the  beginning  of  the 
following  path  to  the  end  of  previous  one.  Such  arrays  were  formed  for  all  isotopes 
and  multifractal  spectra  were  computed.  The  results  of  this  analysis  are  shown  in 
Fig.  6.  The  comparison  of  results  of  both  methods  applied,  to  Irtysh  data,  points 
at  good  concurrence  of  isotope  discrimination  on  their  scaling.  The  qualitative 
behaviour  f{a)  and  7(1/)  (Fig.  7)  curves  are  the  same.  Man-made  isotope  Cs  has 
demonstrated  the  most  large  Holder  exponent  a ~ 1.7  in  multifractal  spectra.  It 
means  that  there  are  hot  spots  of  Cs  contamination  which  form  a sparse  set  having 
dimension  f{a)  < 0.2. 

4 Conclusion 

All  the  results  of  the  multifractal  analysis  that  we  obtained  support  the  existence 
of  multifractal  nature  in  the  terrestrial  radionuclides  contamination  in  the  two 
investigated  regions  and  Irtysh  area.  The  analysis  made  it  possible  to  distinguish 
radionuclide  isotopes  K,  Th,  U,  Cs,  specially  for  measurements  made  on  Irtysh  test 
site. 

Using  the  method  of  computation  topology  for  the  case  of  measurements  paving 
the  territory,  a pre-multifractal  spectrum  of  dimensions  (indexes  of  disconnected- 
ness) for  Irtysh  area  was  calculated.  It  was  found  that  the  behavior  of  disconnect- 
edness index  of  natural  isotopes  differs  strongly  from  Cs  one.  Contamination  of 
natural  radionuclides  exhibits  7 2 for  low  values  and  dimension  decreases  quickly 

with  increasing  of  level.  Spectrum  of  dimensions  for  Cs  points  out  the  existence  of 
’’hot  spots”  of  contamination. 

Taking  into  account  the  multifractal  properties  of  man-made  Cs  found  for  Irtysh 
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area,  it  can  be  pointed  out  the  existence  of  dangerous  for  human  health  anomalies 
on  all  contaminated  territories.  These  ’’hot  spots”  cannot  be  detected  by  using 
spectrometrical  measurements  from  air  accepted  in  Kazakhsatan.  The  results  ob- 
tained can  be  the  background  of  new  methods  of  diagnostics  of  large  territories. 
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